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A reduction of the mass of the rf (958) meson may signal restoration of the Ua(X) 
symmetry in a hot and dense hadronic matter, corresponding to the return of the 
9th, "prodigal" Goldstone boson. We report on an analysis of a combined PHENIX 
and STAR data set on the intercept parameter of the two-pion Bose-Einstein corre- 
lation functions, as measuremed in y/s NN = 200 GeV Au+Au collisions at RHIC. 
To describe this combined PHENIX and STAR dataset, an in- medium n' mass 
reduction of at least 200 MeV is needed, at the 99.9 % confidence level in a broad 
model class of resonance abundances. 



1. Introduction 

Although the quark model exhibits a U(3) chiral symmetry in the limit of 
massless up, down and strange quarks, and in principle 9 massless Gold- 
stone modes are expected to appear when this symmetry is broken, only 
8 light pseudoscalar mesons are observed experimentally. This puzzling 
mystery is resolved by the Adler-Bell-Jackiw Ua{^) anomaly: instantons 
tunneling between topologically different QCD vacuum states explicitely 
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break the £/a(1) part of the U(3) symmetry. Thus the 9th Goldstone boson 
is expected to be massive, and is associated with the 77' meson, which has 
a mass of 958 MeV, approximately twice that of the other pseudoscalar 
mesons. 

In high energy heavy ion collisions at RHIC, a hot and dense medium 
is created. Recent measurements of the direct photon spectrum in y / s^" 
= 200 GcV Au+Au collisions indicate 1 , that the initial temperature in 
these reactions is at least 300 MeV, while hadrons as we know them may 
not exist above the Hagedorn temperature of Tu s=s 170 MeV 2 . Thus 
the matter created in heavy ion collisions at RHIC is hot enough to be a 
quark-gluon plasma 1 . Detailed analysis of the properties of this matter 
indicate that it flows like a perfect fluid 3 , and scaling properties of the 
elliptic flow indicate scaling with the number of constituent quarks 4 , hence 
this matter is sometimes referred to as a strongly interacting Quark-Gluon 
Plasma (sQGP) 5 , or, in more direct terms, a perfect fluid of quarks 3 . 

After this perfect fluid of quarks rehadronizes, a hot and dense hadronic 
matter may be created, where the C/a(1) symmetry of the strong interac- 
tions may temporarily be restored 6 ' 7 ' 8 . Recent lattice QCD calculations 
indicate that such chirally symmetric but hadronic matter may exist below 
the critical temperature for quark deconfincment 9 . In such a medium, the 
mass of the i]' (958) mesons may be reduced to its quark model value of 
about 500 MeV, corresponding to the return of the "prodigal" 9th Gold- 
stone boson 7 . Here we report on an indirect observation of such an in- 
medium 77' mass modification based on a detailed analysis of PHENIX and 
STAR charged pion Bose-Einstein correlation (BEC) data 1041 . 

The abundance of the 77' mesons with reduced mass may be increased 
at low pt, by more than a factor of 10. One should emphasize that the 77' 
(and 77) mesons almost always decay after the surrounding hadronic matter 
has frozen out, due to their small annihilation and scattering cross sec- 
tions, and their decay times that are much longer than the characteristic 
5-10 fm/c decoupling times of the fireball created in high energy heavy ion 
collisions. Therefore one cannot expect a direct observation of the mass 
shift of the 77' (or 77) mesons: all detection possibilities of their in-medium 
mass modification have to rely on their enhanced production. 

An enhancement of low transverse momentum 77' mesons contributes 
to an enhanced production of soft charged pions mainly through the 
77' — > 7] + 7r + + 7r~ — > (ir + + 7T° + 7r~) + 7r + + 7r~ decay chain and also 
through other, less prominent channels. As the 77' decays far away from 
the fireball, the enhanced production of pions in the corresponding halo 
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region will reduce the strength of the Bose-Einstein correlation between 
soft charged pions. The transverse mass (tot = \/m? + p\) dependence of 
the extrapolated intercept parameter A» of the charged pion Bose-Einstein 
correlations was shown to be an observable that is sensitive to such an 
enhanced i]' multiplicity, as pointed our first in Ref. 14 and discussed in 
Ref. 15 . The predicted decrease of A*(tot) data at low transverse mass has 
been observed both by PHENIX 10 and STAR n < 12 at RHIC. 
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Figure 1. The energy and system size dependence of the relative intercept parameter 
in the NA44 S+Pb and STAR Cu+Cu and Au+Au datasets. 



Figure 1 indicates that a "hole" at low transverse mass, in the region 
characteristic to the rj — » rj + n + + n~ — > r/ — > (n + n n~) + 7r + + n~ decay 
chain is not present in the S+Pb data set at CERN SPS at y/s NN = 19.4 
GeV energy 13 , however, it is present and becomes slightly deeper as the 
system size and the colliding energy is increased from Cu+Cu to Au+Au 
and from ^/s^ = 62 A GeV to 200 GeV 12 . 
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2. Modeling and analysis method 

Our main analysis tool was a Monte-Carlo simulation of the transverse 
mass dependence of the long lived resonance multiplicities including the 
possibility of an enhanced rj' production at low transverse momentum, due 
to a partial in- medium Ua{^) restoration and a related rj mass modification. 
This model and the related reduction of the effective intercept parameter of 
the two-pion Bose-Einstein correlation function was proposed first in ref. 14 
and detailed recently in Refs. 16 > 17 > 18 . 

In thermal models, the production cross sections of the light mesons are 
exponentially suppressed by the mass. Hence one expects about two orders 
of magnitude less rj mesons from the freeze-out than pions. This suppres- 
sion, however, may be moderated as a consequence of a possible rj mass 
reduction, and the 77' mesons may show up in an enhanced number. The 
number of in-medium 77' mesons is calculated with an improved Hagcdorn 
formula yielding the following 77' enhancement factor: 

m„i \ "V "V 

This formula includes a prefactor with an expansion dynamics dependent 
exponent a s=s 1 — d/2 for an expansion in d effective dimensions 19 . As a 
default value, a = was taken 14 and, for the systematic investigations, 
this parameter was varied between —0.5 < a < 0.5. Other model param- 
eters and their investigated ranges are described as follows: T conc i in the 
above formula corresponds to the temperature of the medium when the 
in-medium modified rj mesons are formed; its default value was taken to 
be T conc i = 177 MeV 14 and varied systematically between 140 and 220 
MeV. Resonances with different masses were simulated with a mass depen- 
dent slope parameter T e ff — T FO + m(u T ) 2 , where the default values of 
Tfo = 177 MeV and (ut) — 0.48 20 were utilized and systematically varied 
in the range of 100 MeV < T FO < 177 MeV and 0.40 < (u T ) < 0.60 . 

Once produced, the v( is expected to be decoupled from other hadronic 
matter, since its annihilation and scattering cross sections are very small 7 . 
If the T)' mass is reduced in the medium, the observed rj' spectrum will 
consist of two components. If the pr of the 77' is large enough, it can get on- 
shcll and escape. This will produce a thermal component of the spectrum. 
Energy conservation at mid-rapidity implies m*, 2 +Pt, v > 2 = m v > 2 +pT, n ' 2 - 
(In the latter equation the quantities marked with an asterisk denote the 
properties of the in- medium 77', while the ones without an asterisk refer to 



the free 77'.) On the other hand, r/-s with p^ v , < \/m*, — m^ 2 will not 
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be able to leave the hot and dense region through thermal fluctuation since 
they cannot compensate for the missing mass 7 ' 8 , and thus will be trapped 
in the hot and dense region until it disappears. As the energy density of 
the medium is dissolved, the effect of QCD instantons increases and the 
trapped rj' mesons regain their free mass and appear at low px- 

In our recent works of Refs. i 6 * 17 * 1 ^ W e improved on earlier simulations 
of Rcf. 14 , that considered the trapped i]' mesons to leave the dissolving 
medium with a negligible pt- That earlier approach resulted in a steep hole 
in the extrapolated intercept parameter A*(tot) at a characteristic trans- 
verse mass of ttit < 250 MeV 14 > 15 ' 21 . In that simplified scenario the only 
free parameter was the in-medium rf mass, determining the depth of the 
observed hole. In a recent analysis, summarized here, the rj'-s from the 
decaying condensate were given a random transverse momentum, following 
Maxwell-Boltzmann statistics with an inverse slope parameter B^ 1 , which 
was necessary to obtain a quality description of the width and the slope 
of the A*(m T ) data of PHENIX and STAR in the m T w 300 MeV region. 
Physically, B^ 1 is limited by Tpo, so the trapped rj'-s may gain only mod- 
erate transverse momenta. Hence, the enhancement mostly appears at low 
Pt 6 ' 7 ' 8 just as in the first simulations. However, now the slope of "hole" 
of the A*(mx) curve is determined by B" 1 , and, for certain values of the 
model parameters, the data can be reproduced quantitatively. (The A* 
values, actually used in the presented analysis, and their total errors are 
discussed in details in Ref. 18 . Here A™ ax is the A,(wt) value taken at 
hit = 0.7 GeV, with the exception of the STAR data, where the data point 
at the highest hit = 0.55 GeV is considered. Note that the tot dependency 
of the A*(tot) measurements in the 0.5-0.7 GeV region is very weak.) 

We have investigated a broad class of models of resonance production, 
including two different models that produce resonances without assuming 
local thermalization: FRITIOF 22 and UrQMD 23 . Resonance decays, in- 
cluding decay chains, were simulated with JETSET 7.4 24 . 

The FRITIOF 22 Monte Carlo model, based on superposition of nucleon- 
nucleon collisions and the Lund string fragmentation model, cannot de- 
scribe the behavior seen in \<,{m^) / \™ ax even when an arbitrary mass 
modification is considered. On the other hand, hadronic cascade based 
UrQMD 23 , as well as the quark coalescence model ALCOR 25 and the 
thermal resonance production models of refs. 26 > 28 > 27 ; provide a successful 
fit in a certain range of the in-medium rf masses. The main difference 
between the thermal models that we utilized was in those resonance multi- 
plicities that are not yet measured well: ref. 26 predicts a factor of 1.6 more 
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Figure 2. The transverse mass dependence of the relative intercept parameter in the 
PHENIX and STAR datasct is reproduced with an in-medium mass modification of the 
rj mesons using two different resonance models as input. The same resonance models, 
but without in-medium mass modification, cannot explain these datasets. 



rj-s and a factor of 3 more rj'-s than the models of ref. 27 - 28 . The relevant 
resonance fractions of these models are detailed in Table V of ref. 18 . 

The dotted line in Fig. 2 indicates a scenario without an in-medium rj' 
mass reduction, while the dot-dashed and solid lines show the enhancement 
required to describe the dip in the low wit region of A* corresponding to 
the resonance multiplicities of Refs. 26 - 27 , respectively. 

Based on extensive Monte-Carlo simulations, x 2 of the fits to the data 
of Fig. 2 was computed as a function of m*, and B^ 1 for each resonance 
model and each fixed value of model parameters of a, T con d, Tpo and 
(ut)- The best values for the in-medium mass of rj' mesons are in, or 
slightly below, the range yj 'i(2m|- + to 2 ) < m*, < y / 2m'^ < - — m 2 predicted 
in ref. 7 , while all are above the lower limit of m*, > ^/Zvn^ given by 
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Figure 3. Standard deviation contours on the (B~ 1 , m* ,) plain, obtained from 
AtfmxJ/A™" 1 of Monte Carlo simulations based on particle multiplicities using two 
different models for hadronic resonances. The region between the horizontal solid lines 
indicates the theoretical range predicted by quark model considerations, while the dotted 
horizontal line stands for Weinberg's lower limit. 

ref. 29 . The A*(mx) / 'A™ aa: simulations for the best fits of two characteristic 
models are compared to the no-mass-drop scenario on Fig. 2, while the 
1, 2 and 3-cr parameter boundaries are indicated in Fig. 3. Those models 
that describe both PHENIX and STAR A* (m T )/\™ ax data in a statistically 
acceptable manner with the assumption of a sufficiently large in-medium 
T]' mass reduction are all used for the estimation of systematics. The key 
parameters of the best fits are listed in Table 1. 

3. Results 

We have used different input models and setups to map the parameter 
space for a twofold goal: to determine, at least how big rj' in-medium 
mass reduction is needed to be able to describe these datasets, and also to 



I 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


i i i i i i i i i i i - 


= Contours a=0 










Kaneta et al. | 


Kapusta et. al. range : 


z Letessier et al. 

- i i i i l i i i i l i i i i l i i 


Weinberg limit : 

i I i i i i I i i i i I i i - 


) 50 100 150 


200 250 300 



December 23, 2010 1:20 Proceedings Trim Size: 9in x 6in 



cs-gr80 



8 

determine, what are the best values of the in-medium mass modification of 
the rf mesons. Utilizing our indirect method, we have also reconstracted 
the transverse mass dependent spectrum of these rj mesons. 

3.1. Lower limit on the in-medium rj' mass reduction 

We excluded certain regions where a statistically acceptable fit to the data 
is not achievable, thus we can give a lower limit on the rj mass modification. 
At the 99.9 % confidence level, corresponding to a more than 5-ct effect, at 
least 200 MeV in-medium decrease of the mass of the 77' (958) meson was 
needed to describe both STAR 0-5 % central and PHENIX 0-30% central 
Au+Au data on A*(m T ) / \T ax in t/snn = 200 GeV Au+Au collisions at 
RHIC, in the considered model class. 

3.2. Best value of the in-medium n' mass reduction 

We have determined the best values and errors of the fitted m*, and B~ x 
parameters. The best simultaneous description of PHENIX 10 and STAR 11 
relative intercept parameter data is achieved with an rj mass that is dra- 
matically reduced in the medium created in central Au+Au collisions at 
RHIC from its vacuum value of 958 MeV to 340 +g° ± 45 MeV. The 
first error here is the statistical one determined by the 1-er boundaries of 
the fit. The second error is from the choice of the resonance model and 
the parameters (a, T con di Tpo and (ut)) of the simulation. The third er- 
ror is the systematics resulting from slightly different PHENIX and STAR 
centrality ranges, particle identification and acceptance cuts. These effects 
have been estimated with Monte-Carlo simulations, detailed in ref. 18 , not 
to exceed 9.8%, 7% and 3% respectively. The main source of systematic 
errors is the choice of the resonance models. This is due to the unknown 
initial rj multiplicity, hence models like ref. 26 with larger initial 77' abun- 
dances require smaller in-medium rj mass modification, as compared to the 
models of ref. 28 < 27 . 

3.3. Transverse mass spectrum of n' mesons 

In addition to the characterization of the in-medium rj mass modification, 
the transverse momentum spectra of the r\ and rj mesons have also been re- 
ported in Ref. 16 . Fig. 4 indicates the reconstructed spectrum of rj mesons 
m y/ s NN = 200 GeV Au+Au collisions, for simulations based on resonance 
abundances of Refs. 27 > 26 . Normalization was carried out with respect to 



December 23, 2010 1:20 Proceedings Trim Size: 9in x 6in 



cs-gr80 



the rf multiplicity of the model described in Ref. 26 . The spectrum of 
Fig. 4 features a characteristic low transverse momentum enhancement. 
Although PHENIX measured before the r\ spectrum in the pr > 2 GeV 
region 30 , as far as we know the spectrum of the rf particles has not been 
determined before in ^/s^T = 200 GeV Au+Au collisions at RHIC. 

The restoration of the Ua(1) symmetry and the symmetry between mass 
of the T] and the rf mesons is illustrated on Fig. 5. 
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Figure 4. The transverse mass dependent spectrum of the r/ mesons, obtained using two 
different resonance models as input. The band indicates the systematic error, obtained 
from varying the resonance models as discussed in the text. 



4. Discussion 

Detailed analysis of the STAR and PHENIX A* (m T )/\™ ax dataset recorded 
at 7.7, 9.2, 11.5, 39 and 62.4 GeV during 2010 has just been started 31 , 
marking the beginning of the RHIC energy scan program. 
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Figure 5. The left hand side plot indicates the 9 pseudoscalar mesons, the size of the 
peebles is proportional to the mass, respectively. The central plot indicates that the rf 
meson reduces its mass in a hot and dense hadronic medium created in Au+Au collisions 
at RHIC. The right-hand side plot indicates that after the Ua(1) symmetry is effectively 
restored, the mass of the rj and the if mesons will be similar. 



Table 1. Best fits of m*, and B 1 for different resonance multiplicity models, 
followed by x 2 /iVD F and the corresponding confidence level (CL). The integrated 
rj' and rj enhancement factors /„/ and fn are followed by the 5-ct limits of maximum 
in-medium masses. Errors on m*, values represent l-cr boundaries, while the 5-cr 
limits include systematic errors too. 



Resonance m* , B 1 \ 2 „, f„i /17 5-cr limit 

model (MeV) (MeV) NDF ^ m* , (MeV) 



ALCOR 25 


480 ±«5 
530 ±» 


12 


20.2/11 (4.29) 


43.4 


5.25 


< 700 


Kaneta 26 


55 


22.8/11 (4.12) 


25.6 


3.48 


< 730 


Letessier 27 


S40 +50 


86 


18.9/11 (6.35) 


67.6 


4.75 


< 570 


Stachel 28 


S40 +50 


86 


18.9/11 (6.38) 


67.6 


4.97 


< 570 


UrQMD 23 


400 +50 


86 


18.9/11 (6.14) 


45.0 


7.49 


< 660 



At present, detailed data are available from the NA44 collaboration at 
y/s NN — 19.4 GeV 13 as well as from the STAR collaboration at y/s NN = 
62.4 and 200 GeV Cu+Cu and Au+Au collisions, the latter at different 
centrality classes within the 0%-80% range 12 . The NA44 data at ^s NN = 
19 A GeV does not feature an rj' mass drop effect. A positive sign of the rj' 
mass modification is apparent in each case of the STAR datasets, indicating 
that the mass modification effect is nearly at maximum in y/s NN = 200 GeV 
Au+Au collisions and reduces with decreasing centrality, colliding energy 
and system size. We have estimated the magnitude of the system size 
and energy dependence between 62.4 GeV Cu+Cu and 200 GeV Au+Au 
collisions to be not larger than 15%, which is substantially less than the 
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dominant systematic error coming from the choice of the resonance model. 

The dilepton spectrum has been measured recently in minimum bias 
Au+Au collisions at y/s NN = 200 GeV, and a large enhancement was ob- 
served in the low invariant mass region m cc < 1 GeV 32 . Low transverse 
mass enhancement of the r/' and n production results in dilepton enhance- 
ment just in this kinematic range 7 . Estimations using the enhancement 
factors in Table 1 indicate that the observed in-medium rf mass drop is 
indeed a promising candidate to explain this dilepton excess. 

PHENIX recently reported a two-component transverse momentum 
spectrum in dilepton channel direct photon measurements 32 , which provides 
an additional testing possibility to constrain the two component structure 
of the v' spectra reported here. 



5. Summary 

Our report presents a statistically significant, indirect observation of an in- 
medium mass modification of the rj mesons in ^Js NN = 200 GeV Au+Au 
collisions at RHIC. These results were recently published in Refs. 16 > 17 > 18 . A 
similar search for in-medium vj mass modification provided negative result 
in S+Pb reactions at CERN SPS energies 14 . More detailed studies of 
the excitation function, the centrality and system size dependence of the 
A„(m T )/A™ aa: could provide important additional information about the 
onset and saturation of the partial U a{X) symmetry restoration in hot and 
dense hadronic matter. Studies of the low-mass dilepton spectrum and 
measurements of other decay channels of the rj meson may shed more light 
on the reported magnitude of the low pt rf enhancement and the related 
Ua{1) symmetry restoration in high energy heavy ion collisions. 
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